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We derive a general effective many-body theory for bosonic polar molecules in strong interaction 
regime, which cannot be correctly described by previous theories within the first Born approxima- 
tion. The effective Hamiltonian has additional interaction terms, which surprisingly reduces the 
anisotropic features of dipolar interaction near the shape resonance regime. In the 2D system with 
dipole moment perpendicular to the plane, we find that the phonon dispersion scales as yjpf in the 
low momentum (p) limit, showing the same low energy properties as a 2D charged Bose gas with 
Coulomb (1/r) interactions. 

I. INTRODUCTION 

Recent developments in the trapping and cooling of chromium atoms [l[ and polar molecules Q provide a new 
direction for investigating quantum states resulting from the anisotropic dipole interaction. Dipolar effects on the 
condensate profile [3j and elementary excitations [J] have been extensively studied both theoretically and experi- 
mentally. Several exotic many-body states resulted from dipolar interactions are also proposed However, most 
theoretical works so far are based on the pseudo-potential developed by Yi and You [12] within the first Born approx- 
imation (FBA). As a result, these results become not justified when applied to the systems of polar molecules, which 
can have large electric dipole moments and hence strong dipolar interaction to renormalize the scattering amplitude 
beyond the FBA. 

In Ref. 0, Derevianko extended Huang and Yang's approach Q to the anisotropic dipolar interaction, and for 
the first time shew how it may be possible to go beyond the Born approximation in a dipolar gas system. The 
derived pseudo-potential, however, is non-hermitian in the low energy limit, and therefore cannot be easily used for 
constructing the effective theory of strongly interacting dipolar gases. In Ref. [8], the authors studied the systems 
of bosonic dipoles via Monte Carlo calculation, and found that the ground state energy can be well-explained by the 
Yi and You's pseudo-potential within FBA if using a dipole-dependent s-wave scattering. Results of Ref. [8| can 
certainly be applied to a regime of stronger dipole momentum (beyond the valid regime of Yi and You's original 
pseudo-potential), because the higher order renormalization of the s-wave scattering amplitude has been included. 
But its validity to apply to polar molecules with large dipole moment in strong field is still questionable, because 
the higher order renormalization to the scattering amplitude (FBA is the first order perturbation) of the non-s-wave 
scattering channels are not included at all. As a example, in Ref. 0, Deb and You found that the scattering matrix 
element between s-wave and d-wave channels also has strong deviation from their weak interaction result when near 
the shape resonance. It is reasonable to believe that there will be such kind of deviation from the FBA results in 
other channels in stronger dipole moment regime, as usually considered in the polar molecule systems. Therefore 
developing a correct and widely-applicable pseudo-potential and the associated many-body theory for systems of 
dipolar atoms/molecules is still a very important and crucial step for future theoretical and experimental studies. 

In this paper, we derive a complete effective many-body theory which can correctly describe bosonic polar molecules 
both in the weak and strong interaction regime and/or near the SR. The resulting effective Hamiltonian is modified 
by additional three-point and four-point interactions, which may significantly change the nature of condensate pro- 
file/dynamics. For example, when the dipole moment is near the first s-wave shape resonance regime, we find that 
the additional interaction (new terms beyond the FBA) can reduce the anisotropy of the condensate profile. In a 2D 
uniform system with the dipole moment perpendicular to the plane, we find that the phonon dispersion scales as yjqj 
(instead of |q| in the typical Bogoliubov mode) in the long wavelength limit, showing the same low energy physics as 
2D charged bosons with Coulomb (1/r) interactions [lOj. As a result, our theory is important not only in the study of 
strongly interacting polar molecules, but also in the possible application of simulating the liquid phase of 2D charged 
bosons by neutral particles. Such simulation cannot be done in ion traps because of the strong Coulomb potential 
compared to the kinetic energy. These results may be useful in studying the properties of High T c superconducting 
thin film, where the coherent length of Cooper pairs are known to be very small as a composite charged boson [ll| . 

The paper is organized as following: In Sec. HH we first discuss the general scattering theory of dipolar interaction 
and briefly review theories used in previous work. In Sec. IHII we used a exactly solvable model to exam the validity of 
the first Born approximation of dipole interaction. In Sec. IIV1 we derive the correct effective theory and the associated 
Gross-Pitaivskii equation for bosonic polar molecules in 3D system. We then discuss the condensate wavefunction by 
using variational method in Sec. [V] Finally, we extend the 3D results to develop an effective theory in quasi-2D system 
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in Sec. IVI1 and calculate the phonon mode dispersion as well as the Kosterlitz-Thouless transition temperature. We 
then summarize our work in Sec. IVIII 

II. LOW ENERGY SCATTERING THEORY OF DIPOLES 

For the convenience of later discussions, we first briefly review the recent progress on the scattering problem of 
dipolar gas, where the electric/magnetic dipole moment is polarized by the external electric/magnetic field along z 
direction. The most general form of the scattering amplitude between two identical particles in such situation can be 
expressed to be 

/(k,k') = 4vr^^C(fc)r^(fc)y i , m ,(fc'), (i) 

Irn I'm' 

where the scattering matrix element, t\™ (k), depends on the relative incident momentum, k, and the summation 
is over even I for bosons and odd I for fermions. Yi m (f) is the spherical harmonic function of unit vector f. At 
large distances, the inter-particle potential is dominated by the dipolar term, Vd(r) = D 2 (l — 3 cos 2 #)/|r| 3 , where 
D is the electric dipole moment in c.g.s. unit (for simplicity, here we use electric dipoles to formulate the theory 
for polar molecules, while an similar version for magnetic dipolar gas can be also obtained easily); 6 is the angle 
between the distance r and the dipole direction (polarized in z direction). However, at short distances the potential 
becomes much more complicated due to the Coulomb and spin exchange interaction between electrons. Deb and You 
@ first calculated the cross section within a certain model potential and studied how they are changed near the shape 
resonance regime. Based on the numerical results, Yi and You [12j then proposed a pseudo-potential: 

V ps (r) = ^^S(r)+V d (r), (2) 

to calculate the low energy scattering matrix element within the first Born approximation away from the shape 
resonance regime. Here a s — — ^oo(O) ^ s the s-wave scattering length in zero field limit. Within the FBA [Hi]], 
/s(k, k') = J dr e J ( k ~ k ')' r Vp S (r) = — a s — ^a d P2 (cos#k-k')j where a d = MD 2 /h 2 is a length scale and P%(x) is 
the Legendre Polynomial. #k is the angle between the momentum k and z axis. As a result, the associated matrix 
elements become: t B ^™' (k) = 4^a d i l ~ 1 ' / dQ, r Yi rn (f)Y l 1 m ,{f)Y 20 {f) ^ji(r)ji>(r) with ji(r) being the spherical 
Bessel function. Here r c is a cut-off in the atomic length scale and therefore we can always take kr c — > in the 
low energy limit. In above FBA result, all the short-ranged effects are included in the s-wave part (a s ) only, while 
all other matrix element, t_B/ r ™ (0), are proportional to the same length scale, ad oc D 2 . Therefore, it is easy to see 
why such results cannot be valid when the dipole moment (or external field) is sufficiently strong as higher order 
renormalization becomes important. 

In Ref. [6(, Derevianko developed a different pseudo-potential for dipolar interaction to go beyond the FBA. 
Although the most general expression of the pseudo-potential is derived for each scattering channels and the results 
are in principle applicable to strong interaction regime, but only one terms (the scattering between s-wave and d-wave, 
i.e. the too(0) = *2o(0) term in Eq. {!])) is evaluated within the leading order perturbative method (equivalent to the 
FBA level). In fact, we observe that Derevianko's result for the on-shell scattering channel (|k| = |k'|) is equivalent 
to the first two terms of the FBA result (i.e. the full /s(k, k') is replaced by — a s — ^ [^(cos^k) + ^(cos^k')]) 
using Yzo(k) = ^/5/47rP2(cos 9k))- This explains why the meanfield calculation by Yi and You (Ref. [13], which 
included only the s-wave and s-d scattering channel of the pseudo-potential of Ref. (§]) is not consistent with the 
previous result even in the weak dipole moment regime, where the FBA is supposed to be valid. (We note that this 
inconsistence still exist even if Yi and You have ever used the corrected coefficient derived by the Erratum of Ref. 
Q. The key point is that contributions from all other scatteringchannels are all proportional to D 2 within the FBA 
and hence cannot be neglected compared to £qo (0)0 ^ n Bortolotti et al. claimed that V ps (r) in Eq. @ 

can be a good pseudo-potential if only one uses a dipole-dependent s-wave scattering length (i.e. a s (D)). However, 
their results cannot apply to the strong dipole moment regime when the shape resonance occur in other (different 
from s-wave) scattering channels due to the complicated electronic density distribution and/or spin exchange effect 
in a realistic polar molecule. Therefore, a general and useful approach to study the low energy many-body physics of 
strongly interacting polar molecules is still needed. 

III. CRITERION FOR THE FIRST BORN APPROXIMATION 

For completeness, now we explicitly examine the criterion for justifying the FBA in the low energy limit. We 
consider the following model potential: Vmdi(r) = Vd(r) for |r| > r c , and Vmdi(r) = oo for |r| < r c . Although 
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FIG. 1: Series expansion for effective interaction in the ladder approximation. Solid line represents Green's function of bosonic 
particles and zig-zag line is for bare interaction. Here p = |(pi — P2) and p' = ^(ps — P4) are the half the relative momentum, 
and P = pi + P2 = P3 + P4 are the total momentum of the two scattering particles with frequency Pg (see also Ref. Hi 



this model potential is over-simplified compared to the realistic interaction potential between polar molecules, it still 
catches the most important feature, anisotropic dipolar interaction, and hence should be useful in studying the validity 
of the Born approximation in the low energy limit. The full scattering wavefunction if)(r), can be solved by: 
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where j dr' is for |r'| > r c only, and 



V-o(r) = 4vr^ 



cos 5i (k)jni (k, r)Y{^ (k) Y im (f ) 



(3) 



(4) 



is the exact scattered wavefunction for the hard core potential of radius r c without dipole moment. Here we have defined 
jni{k,r) = ji{kr) — tan Si(k)ni(kr) with Si(k) = tan" 1 (ji(kr c ) /ni(kr c )) being the scattering phase shift. ji{x) and 
ni(x) are the conventional spherical Bessel functions. The Green's function, G(r, r'), satisfies V 2 G(r, r') + k 2 G(r, r') = 
— 47r^(r — r') with the boundary condition G(r c r,r') — 0, and therefore can be evaluated by using separation of 
variables. After some straightforward calculation, the Green's function can be expressed to be 



G(r,r') = J2 Y ? m ( 
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where ?">(<) is the larger(smaller) one of r and r' , and h\ ; (x) = ji(x) + ini(x). 

Within the FBA, the scattered wavefunction is given by the first order iteration: V's(r) = ?Ao( r ) — 
W /' T~^( t j r ') 1 ' / mdi( r ')V'o( r ')- Therefore its validity relies on the assumption that the change of the wavefunction is 
much smaller than ipo{ r ) i n the whole range of space [14| . We can therefore define a parameter, £, to measure the 
deviation of V'B( r ) : £ = limt^o lim| r |^ rc |i/ , s( r ) — "00 ( r ) | / 1 "00 ( r ) I - I n such limit, we have i^o( r ) ~ + 0(kr c ), where 
Ar = |r| — r c -C r c . Expanding the Green's function, G(r, r'), in the small r regime, we obtain ipsi^) — "0o( r ) = 
T 1 ' f^f ' ^77 + ( fcr c)- As a result, the condition to justify the FBA is £ = f^f^ < 1 [Hi- For example, we 
consider the magnetic dipolar atom, 52 Cr, with r c ~ lOOao as the typical length scale of van der Waals interaction. 
We find £cr ~ 0.4 < 1 and this explains why results obtained in the FBA for 52 Cr are comparable to experiments 
3.4]. However, for polar molecules with electric dipole moment of the order of a few Debye, the value of £ can easily 
be several hundred or more, where a shape resonance can occur in different channels to breakdown the FBA result. 
Therefore, in order to correctly describe the effective many-body physics of polar molecules, one needs a self-consistent 
theory beyond the pseudo-potential, Vp S (r), and the first Born approximation. 



IV. EFFECTIVE HAMILTONIAN IN 3D SPACE 



To study the low energy physics of a general dipole interaction in the many-body medium, one has to use an 
effective two-particle interaction, T, which is just the vertex function integrating out all the contribution of virtual 
scattering in high energy limit [1 61 ] . A full calculation of the vertex function is usually not available (except in some 
special models of ID systems), but can be well-approximated by using the standard ladder approximation (see Fig. 
[T]). It is well-known that such ladder approximation is correct in the low density limit, and is therefore a very suitable 
approximation for systems of dilute cold atoms/molecules. Following the standard approach to evaluate the Bethe- 
Salpeter equation of bosonic particles [16j . we can calculate the effective two-particle interaction (i.e. vertex function) 
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within the ladder approximation by using the two-particle scattering amplitude, /(p,p') : 

J^rfe, p', p, p ) = -/(p, p') + ^E / Cp, k)/( P ', k)- ( £ + 2W ,L fc2 + t0+ + fc2 _,L z0+ ) - ( 6 ) 

where e = $ (hP - H 2 P 2 /AM) is the total kinetic energy in the center-of-mass frame; \i is the chemical potential 
and Q is the system volume. 

Using the fact, /(p',k)* = /(k, p'), the final term of Eq. © can be evaluated explicitly by integrating over the 
solid angle of momentum k in the scattering amplitude, Eq. Furthermore, since the partial wave scattering 

matrix clement, t\™ (fc), is known to be insensitive to the incident momentum, k, in the low energy limit, we can also 
neglect their momentum dependence and replace their value by a constant, t\™ (0). As a result, the last term of Eq. 
© can be calculated to be 



(4tt) 



k 2 dk 

(2tt) 3 \e + 2Mn/h 2 - k 2 + ir) ' k 2 - p' 2 - ir\ 
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t l ^(0)4Z(0)YC m (p)Y l>m {p') 



(7) 



where we have set m — m! due to the rotational symmetry about the polarization axis (z). It is easy to see that the 
real part of the integration cancels out, and the imaginary part proportional to y / 2M/i/h 2 in the limit of low energy 
scattering (\p'\, e — > 0). Therefore, the last term of Eq. (J6]) can be shown to be negligible when comparing with with 
the second term, /(p, p'), in the low density limit, i.e. (\t\™ (0)rig^| <C 1. Here n 3 D is the 3D particle density). As a 



result, in the low energy and dilute limit, one can use r(p,p') = 



M 



/(P)P') as an effective "pseudo-potential" in 



momentum space (there is no dependence on total momentum and energy in such limit and we could omit them in 
r). Note that, different from the FBA used in the literature, we do not have to assume weak bare interaction in above 
derivation (strong interaction may still give small value of scattering matrix element, \ti™ (0)|, in the low energy 
limit, just as in the usual s-wave scattering of cold atoms). Complicated electronic structure and shape resonance 
effects are all included in the full calculation (or experimental measurement) of the matrix elements, t\™ (0) in all 
channels. In the rest of this paper, we will study the general effective theory and possible new many-body physics 
beyond the FBA without directly evaluating the scattering matrix elements. 

Using the derived pseudo-potential (or effective interaction), T(p,p') = /(p, P')> we can write down the 

interacting Hamiltonian in momentum space by using the second quantization formalism: 



H, 



p,p',p 



Lp +p a T !p_p«ip-p'«Ip + p' r (P: P'). 



(8) 



where a p and a' p are field operator for bosonic polar molecules at momentum p. The momentum summation from 
now on is restricted to low momentum regime as implied by the effective interaction, T(p, p'). In order to address the 
effect of pseudo-potential beyond the FBA (Eq. ([2])), we can divide the contribution of pseudo-potential, T(p,p'), 
into three parts: 



r(p, P ') = 



4nh 2 a s 
M 



^(p-p')-^-/a(p,p'), 



(9) 



where the first term is from the known (dipole moment dependent) isotropic s-wave scattering, the second term is the 
usual FBA result for anisotropic dipolar interaction, the third term, /a, is the scattering amplitude deviated from 
the known FBA results. It can be denoted to be 



/ A ( Pl , P2 ) = -47r]T'^E ^YLWfvM 



(10) 



with Aa 
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f B;rrT(0)) being the difference between a full matrix element and its FBA result. Here 
Note that in the limit of a weak external field, we have following orders of 



0{D 2 ), and Aa 



(m) 



magnitudes: a s . ,, ^ , 

r, shown in Eq. © has a very smooth connection with the known FBA results [3.4,8.12] in the limit of small dipole 
moment. From Eqs. © and it is straightforward to write down the full effective Hamiltonian to describe the low 



0(D 4 ). Therefore, it is easy to see that the pseudo-potential, 
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energy many-body physics of polar molecules: 

H = ^(e p - Ai)aJ,a p + tz E a k a P2^ x t(pi - P2) 



Pl:P2 



i E ^p +p1 «1 



Pl,P2,P 



P-Pl 2 P ~P2 2 P + P2 



4nh 2 a s 
M 



Vd(Pi 



P2J 



4nh 2 a ( 
M 



/a(Pi,P2) 



(11) 



where V^ xt (p) is the external trapping potential in momentum space. Note that Eq. (fTT]) has included scattering 
from all channels and is also consistent with the FBA results in the weak dipole limit (|/a| oc 0(D a ) as D — > 0). 
When the external electric field is strong enough, there will be some modification of the scattering amplitude to be 
beyond the results of first Born approximation even in channels different from s-wave, i.e. Aa\™ ] ^ for I, V ^ 0. 
Calculating the magnitude of such modification beyond the FBA has to be based on the first principle calculation 
of two scattering molecules, and is beyond the scope of this work. Our interest in the current paper is to study the 
effective Hamiltonian and the possible many-body physics when /a is known. 

In order to compare with the existing theory of weakly interacting dipoles 3.4.8.12], it is instructive to express Eq. 
(fTTj) in real space. Details of the transformation is shown in Appendix [A"l The final result is 



H 



drip' (r) 

2^h 2 
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M + v; x t(r) 



dR 



^(r) 



^(R+ -)V> f (R- 



£)£A a <,%(R) 



dr I dr'V ps (r - v')^ (r) ft (r') V> (r') ip (r) 
2ft 2 



h.c. 



M ^ 

i.v 



^"E A < 



dR 



1\ 



;mV R)^ m (R) 
(12) 



where has excluded (1,1') = (0,0), (0,2) and (2,0), and we have used the fact that Aa|™' ) 7^ only for \l — l'\ — 
0, 2,4, • • • due to the anisotropic nature of dipole potential, V^(r) oc Y2o(f), and its higher order effect. As shown 
in Appendix [Aj we have defined 0; m (R) = ^jl'' J dr Yl ™j r ^ V>(R + |)^(R — §) as a "pairing" operator in angular 
momentum (I, m) channel with a spatial "wavefunction" Yi m (f)/r 3 . Although such pairing operators do not represent 
true composite particles, but can be used to describe the relative motion of two dipoles before and after scattering: the 
first term of the last line indicates an association-dissociation process between a pair and two dipoles, while the last 
term describe a transition between "pairs" of different angular momentum channels. These two novel interaction terms 
should bring complete new physics in a strongly interacting polar molecules, and is worthy for further investigation 
in the future. 

Starting from the effective Hamiltonian, Eq. (fT2|) , we can also derive the associated Gross-Pitaeviskii type meanfield 
equation for condensate dynamics by using ihdip/dt = [ip,H] and approximating the bosonic field operator, ip(r), to 
be a c- number, 'J'(r). The resulting equation can be written as following form: 
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dt 



-ft 2 V 2 
2m 

2^h 2 
2(47r) 2 ft 2 



fi + V eKt (r) + 



dr' 
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1)H 



*(r) 



2'/2 



-Y l0 (f') 
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(m) 
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47T 
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(Z + l)!! 7T 



2'/ 2 



4tt 



(2tt 



,3^ 



(I' + 1)!! 7T 

2 1 '/ 2 ri, 



(13) 



where ^(r) = (V>(r)) is the condensate wavefunction. Similarly, one can also derive associated Bogoliubov-de Genne 
equations for the elementary excitations. We note that the effective Hamiltonian, Eq. (fTTj) and Eq. (fP2]) . and 
meanfield equation, Eq. (fl~3|) . contain all the effects beyond the simple FBA results, and they will reproduce the 
known FBA results when taking Aa[™' = 0. 
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V. 3D CONDENSATE PROFILE 
A. Gaussian variational wavefunction 



To study the aspect ratio and the stability regime of the condensate profile, it is convenient to use the variational 
approach [17| . Here we use a Gaussian type trial wavefunction, 



*(r) = 



Nexp(-p 2 /2R 2 - z 2 /2R 2 z ) 



tt 3 / 4 R R 



1/2 



(14) 



for the condensate wavefunction in harmonic trapping potential: V^ xt ( r ) = \Mlo 2 (x 2 +y 2 ) + ^Muj 2 z 2 , where loq and uj z 
are the associated trapping frequencies. Here N is the total number of dipoles, and i?o and R z are the Gaussian radii 
of the condensate in the x — y plane and along the z axis respectively. The variational energy can be obtained easily 
from the effective Hamiltonian in the momentum space, Eq. (fTTj) . via replacing at by $t = (ak) = ^ J dr^f(r) e _lk r . 
We therefore obtain 

E{R ,(3) l + 2(3 2 R 2 (2k 2 + (3 2 ) 2N& s 16 A 2 (/3) Nd, d 



2P 2 R 2 
32/37V 



2k 2 



^RlP 



2irRl 



1=2 



(15) 



where Ai{(3) 



_\M a 

(l + (/3 2 -l)x 2 )3/ 2 • P 



R z /Rq and k = ujo/lo z are the condensate and trapping aspect ratios. 



We have also scaled all the length scales (a s , a d , AafP, 
(i.e. -Ro = -Ro/aosc,o etc.), and used i?o = Nh 2 /2ma 2 sc 
The first two terms in the right hand side of Eq. 



and i?o) by the horizontal oscillator length, a OSCj o = Wh/MuJo 
o as the energy scale. 

(fT5T> are from the kinetic and potential energies respectively, 
and the third is from the s-wave scattering channel. The fourth term in from the contribution within the first Born 
approximation and the second line is from the effects beyond the FBA. Again we find that the whole mcanfield energy 
of Eq. (|15[) will become the same as calculated within the First Born approximation by taking Aa|°^ +2 = 0. Using 

the fact that Ao(/3) = (4/3)" 1 , we find that the contribution of the FBA is of the same form as the term with Aa^l 
(both of them are proportional to A 2 (/?)). However, such coincidence is simply due to the special form Gaussian trial 
wavefunction. Using other trial wavefunctions can easily give different aspect ratio dependence of these two effects. 
Besides, we also note that Ai((3 = 1) = for I ^ 0, showing that for a spherically symmetric condensate (j3 = 1), 
only s-wave scattering channels are relevant: scattering in finite angular momentum channels are cancelled out due 
to spherical symmetry of the condensate profile. When the condensate profile is highly anisotropic due to external 
confinement (say in cigar shape, (3 ^ 1 or in pancake shape, j3 1), the effects beyond the FBA will become very 
crucial. 



B. Example: near shape resonance 



For the general form of effective Hamiltonian of Eqs. (fTTj) and (TT2"j) . values of a,,7^ have to be obtained from the 
first principle calculation [10], which is however beyond the scope of this paper. In fact, due to the highly nontrivial 
inter-molecule interaction in short-distance, the low energy scattering matrix element, t\™ (0), can be very different 
from the results of Born approximation in strong dipole regime. Here we consider the simplest case to study the effect 
beyond the FBA: we assume the external electric field is still weak but near the first shape resonance regime, where 
the shape resonance occurs in the s-wave channel so that both £qo(0) an< ^ ^oo(O) are strongly deviated from results in 
the weak interaction limit. Scattering matrix elements in other channels are less affected because of the centrifugal 
potential for I ^= 0. This picture is also consistent to the numerical results shown in Ref. @, where their numerical 
results of t|o(0) is almost unaffected by the shape resonance in the s-wave channel. (But it does not exclude the 
possibility to have significant deviation in other channels in the regime of much stronger dipolar interaction.) 

Under such assumption, we may consider Aa^ ^ 0, and A[ ^ = for all (I, V) ^ (0, 2), (2,0) or (0,0). As a result, 
the variational energy of Eq. (4) becomes (using Aq((3) — 1/4/3): 



E(Ro,0) 1 + 2/3 2 R 2 {2k 2 +P 2 ) 



Eq 
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(16) 
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From above result, we find the contribution of the Aa ( ' 2 term can reduce (since A<2q ' 2 > near the first shape 
resonance, see Ref. @) the effect of anisotropic feature of dipole interaction. Although this result is derived from the 
Gaussian trial wavefunction, such reduction of anisotropy of the condensate wavefunction should be still qualitatively 
correct for the correct condensate profile. 

We note that the ground state energy and the pseudo-potential study has been also discussed in Ref. [8|], where 
they include the dipole dependence in the s-wave scattering length (i.e. a s (D)) and use the FBA results (Eq. ([2])) for 
the dipole interaction near the first few shape resonance. In other words, they did not consider the effect of strong 
deviation of fgo(O) fr° m the FBA, which is a very significant result as shown in Ref. |9j. Therefore, the large value 
of s-wave scattering length near the shape resonance in Ref. *8j may have smeared out the contribution of Aa^. 
If considering an even stronger dipole moment (larger than the value for the first few shape resonance), where the 
scattering amplitudes may deviate from the FBA result in all channels, one has to solve Eqs. (fTTj) - (fT"3)) with finite 

values of Aajy for the correct many-body physics of polar molecules. 



VI. EFFECTIVE HAMILTONIAN AND EXCITATIONS IN 2D: 



In a 2D homogeneous system, we can assume that the wavefunction in the z axis is of Gaussian type: 4>{z) = 
~ z / 2R z ; where R z is the width of such quasi-2D potential layer. After integrating out the degree of freedom 



in z direction (i.e. along the direction of external electric field) of the 3D effective Hamiltonian, Eq. (fTTj) . we obtain 



pi,P2,q 



ik E ^^WU 



bp bp . 

2--P2 2-+P2 



(17) 



Pl,P2,P 



where we define b p and 6p to be the field operator in 2D system with p being the in — plane momentum vector from 
now on. fij_ is the 2D area, and V s = J^)^ S R is the contribution of s-wave scattering. Vs(q) = (~^ £ ) > 

where g(x) = 1 — (3y/n/2)x e x Erfc(a;) with Erfc(a;) being the complementary error function. We also have 



MR Z 



(18) 



i,v m 



to account the contribution beyond FBA, where 

F lm (p) = R z [ d Pz e~ p i R " Y lm (p) = f°° dx 

J J — OG 



^ 21 + 1 {l-\m\)\ 



4tt {l + \m\)\) 



Pi 



x 2 + p 2 R 2 



im<t>p 



(19) 



with (ftp = tan 1 {p y /px) being the angle in 2D plane. 

At zero temperature, the dipolar atoms/molecules condense at p = 0, so that the total energy is E 2 d = 
2H7 Ws + V B (0) + Va(0, 0)], with the chemical potential being u = n 2D [V s + V B (0) + V A (0, 0)], where n 2D = N/Qj_ 

is the particle density in the 2D plane. Keeping only the condensate part (S = = y/n^o) and the quadratic order 
of fluctuations (p ^ 0), the effective Hamiltonian become: 



H, 



eff 



E ( £ - M ) a p°p + 9?T7 E ^ + Vb ^ + V ^ °)1 ( a P a -P + a -P a : 



N 



2f2 



+^ e h + 2v ^ + 2v ^ + 2v - (!< f ) + 2v - (I -f 



(20) 



where we have used the fact that Va(pi,P2) = Va(— Pi,— P2) = Va(P2,Pi)- Finally, we could use the Bogoliubov 
transformation to diagonalize above Hamiltonian and obtain the following phonon excitation spectrum: 



f - + 2n 2D (V s + V B (p) + W + (pj) 
Am 



(21) 



8 




where W±(p) = Va (§, §) + F A (§, -f ) - Va(0, 0) ± Va(p, 0) accounts the effects beyond the FBA results. 

Similar to the 3D case, now we study the situation when only Aajj ^ 7^ in Va(Pi,P2), an d obtain Va(Pi,P2) = 

^tiFL^ 1 "^ Q o°2 [.9(|Pi|^z) + <?(|P2|-Rz)]- The calculated dispersion, u> p , for different values of Aa^l > are shown in 
Fig. [2] There are two significant effects to be noted: First, in the short wavelength regime, the roton minimum, 

predicted || as a feature of dipolar interaction for 2D bosonic polar molecules, becomes weaker as a^l becomes 

stronger. Secondly, in the long wavelength limit, instead of the typical linear dispersion Q, we find uj p — Cy/\p\R z (l + 

0(\p\R z )) with the prefactor C = 2V6(57r 3 ) 1 / 4 y // Aa^ (a s + §a d - 2v / 5Aa^) x . As a result, the phase 

fluctuation becomes much stiffer than predicted in the FBA, showing an enhancement of the condensate/supcrfluid 
density at zero temperature. More precisely, we can calculate the normal fluid density, p n , according to the transverse 
current correlation function [181 ]. The sublinear dispersion of lo p gives p n 

( T ) = v i*M 2 c»S 7 ' which shows a much 

smaller temperature (T) dependence than the result obtained for linear dispersion (p n = 3 ^^mk^ ^ ^ p = cifi.|p| 

[III]). According to Landau's two-fluid model and the universal relation between the 2D superfluid density and the 
Kosterlitz-Thouless transition temperature (T c ), the superfluid transition temperature (T c ) of 2D dipolar system is 
then determined by fcsT c = irti 2 p s (T c )/2M = irh 2 (ri2D — p n (T c ))/2M. At temperature below T c , the single particle 
correlation function has a power-law decay with zero condensate density. These results are also equivalent to a 2D 
charged Bose gas (V(r) — Q 2 /r) [10] with an effective charge, Q = C R Z M j n%T>Y? . Such an interesting equivalence 
implies a possibility to use neutral polar molecules to simulate a 2D charged boson system in liquid phase (not 
doable for ion traps due to the strong Coulomb interaction and large atom mass), which may be important to the 
understanding of the superconducting Cooper pairs in High T c thin film [Tlj . 



VII. SUMMARY 



In summary, we have developed a full effective many-body theory for 3D and 2D dipolar Bose gases beyond the 
simple first Born approximation. One of the significant consequence is that the dipolar interaction effect in the 3D 
condensate can be reduced near the shape resonance regime. For the 2D system (highly anisotropic regime), such 
effect brings a significant change of the low energy excitation spectrum. We believe there should be more interesting 
results for a polar molecule system in strong external field regime, where all scattering channels (besides of the ig§ 
channel) can deviation from the FBA significantly. Our results therefore are especially important for the future 
studying of the many-body properties of strongly interacting polar molecules. 



VIII. ACKNOWLEDGEMENT 



We thank G. Baskaran, A. Derevianko, D.-H. Lee, S. Ronen, L. You, and W.-C. Wu for fruitful discussions. Part 
of this work was done in KITP, Santa Barbara. Our work is supported by NSC Taiwan. 



9 



APPENDIX A: EFFECTIVE HAMILTONIAN IN REAL SPACE 

We note that the single particle part, s-wave scattering, and the FBA part of interaction can be easily transformed 
in to real space as shown in the literature, therefore here we just show the results for the interaction part beyond the 
FBA, i.e. the last term of Eq. (fTT|) . To Fourier transform the effective Hamiltonian, we use a p = f -^Lxpir) e~ 4pr 
and the interaction part beyond FBA (denoted by H'j) becomes: 



Hi 



-Airh 2 1 

M m 

(47r) 2 ft 2 1 



M 2Q 3 



Y «lp+p i a |p-p 1 «Ip-p 2 «ip +P2 /A(Pl ! P2) 



Pl,P2,P 



Y ]T'A4 m V'-^Upi)iWp 2 ) 



pi,P2,P ll'ra 



x J dndr 2 dr 3 dr 4 e^+v^e^-v^e-^-v^e 



M 2ft 2 



Y Y' Aa ^ il '~ lY i r m(Pl)Ylm(p2) 



Pi ,P2 11' m 



J dKdv l dv 2 e^e-^ftCR + y )^(R - y )^(R - y )Vi(R + y ), 



(Al) 



where we have integrated out total momentum, P, and center of mass position before scattering, R' = (r3 + Ti)/2. 
Note that the summation of angular momentum quantum number, ^2w m ', has exclude the pure s-wave scattering 
channel, I = I' = m = 0. 

Now we consider the expansion of plane wave in spherical harmonic functions: 



e* k - = 47r£i , j I (fcr)Y^(f)F Im (fc) 



I III 



or equivalently 



dn ke ^Y; m (k) 



4iri l ji(kr)Y l * m (r) for r ^ 
for r = 



(A2) 



(A3) 



Eq. (|A3[) suggested that it is more convenient to separate terms with zero angular momentum quantum number, 
Aa^ = Aa^, , from others in Eq. (|A1[) before carrying out the momentum integral. (We note that the higher order 
correction of dipolar interaction, Vd(r) oc Y2o(r), can couple s-wave to higher moemntum channels in the strong dipolc 
momentum limit.) As a result, Eq. (|Aip can be rewritten to be 



Hi 
H' n 



H n + H' I2 



y J dRdr 1 dr 2 ^ (R+ ^) V 3t (R _^ ) ^ (R _| ) ^ (R+ ^ ) 

tp2 V2 (^o(pi)*bo(p2) + Wi^ofe)) 



^_iy/2 A „(o) 



(47r) 2 ^ 2 
M 

(4tt) 2 H 2 
+ M 

{Airfh 2 



J dRdr^^R + -^(R - y)V>(R)^(R) x £ 1 ^=b» _ £ e^-^*^) 



/ dRdra^CH^WCR- f)^(R+ y) * £ 1 ' 01 ^5>"™y i0 (p 2 ) 

l L P2 



(A4) 



H' I2 = |rfRdr 1 dr 2 v3t(R+^)V3 t (R-^)^(R-|)^(R+|) 



(A5) 
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where J2w m " m the last line is a summation excluding any terms with I or I' = 0. To get Eq. (j A4|l . we have used 
the fact that f2 _1 e lp r = <5(r) and have integrated out one of the relative coordinate. Now we can integrated out 
the solid angle of momentum variables (pi and P2) by using Eq. (|A3|) for I 7^ 0: 



(A6) 



where A is the momentum cut-off in atomic length scale (~ J"" 1 ), due to the nature of effective Hamiltonian obtained 
by integrating out the high momentum/energy contribution within the ladder approximation of Eq. ([6]). In order to 
regularise it to get a universal expression, we can introduce another high momentum cut-off, e~~ ap , inside the integrand 
and taking a to zero (a — > + ) in the final results. Using the fact that only even angular quantum numbers (I = 2m') 
are relevant for the scattering between bosonic polar molecules, and applying the following identity: 



/>oo 

lim / p 2 j2m'{pr)e~ ap dp 



(2m' + 1)!! 7T 



for mf 7^ 0, we can simplify Eq. (|A6|) further and rewrite the effective Hamiltonian {H' n and Hi2 f ) to be 



(A7) 



M 



dUdr 



M 



dR 



ft(R + )^t (R ) £ Aa<^ (R) + h.c. 



(47r) 2 n 2 

W 



H 'l2 



y'rfRdr 1 dr 2 v3 t (R+y)^ t (R-y)^(R-y)^(R 



47r r* fn) (/ + 1)!!7r 



viV"A« (m 

X 2 Z^ ^ a «' (27r)3"' mV 17 2'/2 r? 

^^"^Aa^/dR^L(R)# m (R) 



47T 

(2^ 



-) 

2 7 

(i' + l)!! 7T 



2''/2 r | 



where we have defined an effective pairing operator: 



(I 4. 1 VI C dr - r a r 

Im (R) = / — ll m (f)^(R + -ty(R - -) 



h.c. 



(A8) 



(A9) 



(A10) 



in angular momentum (l,m) channel to simplify the notation. Therefore, after adding back the known single particle 
Hamiltonian and the FBA results together, we can obtain the final total effective Hamiltonian in real space as shown 

in Eq. mr 
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